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Abstract 

We report muon spin relaxation (//SR) and magnetization measurements, to- 
gether with synthesis and characterization, of the Li-intercalated layered su- 
perconductors Lij^HfNCl and Li^jZrNCl with/without co-intercalation of THF 



(tetrahydrofuran) or PC (propylene carbonate). The 3-dimensional (3-d) su- 

perfluid density ns/m* (superconducting carrier density / effective mass), as 

well as the two dimensional superfluid density ns2d/'m*ah (2-dimensional (2- 

d) area density of superconducting carriers / ab-plane effective mass), have 

been derived from the /xSR results of the magnetic-field penetration depth 

Xab observed with external magnetic field applied perpendicular to the 2-d 

honeycomb layer of HfN / ZrN. In a plot of Tc versus ns2d/'^*ah^ most of 

the results lie close to the linear relationship found for underdoped high-Tc 

cuprate (HTSC) and layered organic BEDT superconductors. In Li^jZrNCl 

without THF intercalation, the superfluid density and Tc for x = 0.17 and 0.4 

do not show much difference, reminiscent of /xSR results for some overdoped 

HTSC systems. Together with the absence of dependence of Tc on average 

interlayer distance among ZrN / HfN layers, these results suggest that the 2- 

d superfluid density ns2d/i^*ab ^ dominant determining factor for Tc in the 

intercalated nitride-chloride systems. We also report //SR and magnetization 

results on depinning of flux vortices, and the magnetization results for the 

upper critical field Hc2 and the penetration depth A. Reasonable agreement 

was obtained between /LtSR and magnetization estimates of A. We discuss the 

two dimensional nature of superconductivity in the nitride-chloride systems 

based on these results. 
78.30.Hv, 78.40.Ha, 72.80.Ga, 75.45.+j 
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Typeset using REVT^ 



I. INTRODUCTION 



Layered superconductors, such as high-Tc cuprates (HTSC) or organic BEDT systems, 
have been a subject of extensive research effort for decades. These systems show rich novel 
phenomena, including superconducting fluctuations, pancake vortex, complicated vortex 
phase diagrams, and interlayer Josephson effects. High-Tc cuprate superconductors have 
been investigated extensively as prototypical layered superconductors. The cuprates have 
a merit that their carrier concentration can be controlled by chemical substitutions and/or 
oxygen contents. On the other hand, it has not been easy to control the interlayer distance in 
cuprates. In general, superconductivity of these layered systems is deeply related to in-plane 
features as well as interlayer couplings. For overall understanding of superconductivity in 
these systems, it would be essential to elucidate interplay between in-plane and interplane 
properties. Despite extensive research effort, however, detailed roles of dimensionality are 
yet to be clarified in these systems. 

Recently, superconductivity was discovered in ZrNCl and HfNCl intercalated with alkah 
atoms (Li, Na, K) [1,2]. Systems based on ZrNCl have superconducting transition temper- 
atures Tc < 15 K, while those based on HfNCl have Tc < 25.5 K. These systems can be 
co-intercalated with organic molecules, such as THF (tetrahydrofuran) or PC (propyrene 
carbonate). The parent compounds ZrNCl and HfNCl are insulators which have a layered 
structure as shown in Fig. 1(a). Zr(Hf)-N honeycomb double layers are sandwiched by 
CI layers, and a composite Cl-(ZrN)-(ZrN)-Cl layers form a stacking unit. Adjacent stack- 
ing units are bonded by weak van der Waals force. Alkali metal atoms and polar organic 
molecules such as THF or PC can be co-intercalated into the van der Waals gap of the 
parent compounds as shown schematically in Fig. 1(b). Intercalated alkali metal atoms are 
supposed to release electrons into Zr(Hf)-N double layers, which makes the system metallic 
and superconducting. On the other hand, intercalated organic molecules expand interlayer 
distance without changing Zr(Hf)-N honeycomb double layers. So we can control two sep- 
arate parameters, carrier concentration and the stacking unit distance, in a single series of 



3 



nitride chlorides with the common superconducting slab. This unique feature could allow 
studies of layered superconductors from a new angle. 

In this paper, we will present synthesis and characterization of a series of intercalated 
ZrNCl and HfNCl samples, together with studies of their superconducting properties using 
muon spin relaxation (/xSR), magnetization, and resistivity measurements. A part of this 
work was presented in a conference [3], where we reported /xSR results in HfNCl-Lio.s-THFo.a, 
showed that Tc and the 2-d superfluid density in this system follow the correlations found 
in cuprates and BEDT systems, and discussed that this feature likely comes from departure 
from BCS condensation, which can be understood in terms of crossover from Bose- Einstein to 
BCS condensation. Subsequently, Tou et al. [4] reported an NMR Knight shift study which 
inferred a rather small density of states at the Fermi level in HfNCl-Li-THF, and discussed 
difficulty in explaining the high transition temperature in terms of the conventional BCS 
theory. Tou et al. [5] also reported a rather high upper critical field Hc2{T — > 0) ~ 100 kG, 
for the field applied perpendicular to the conducting planes, from magnetization and NMR 
measurements. 

Extensive //SR measurements of the magnetic field penetration depth A have been 
performed to date in various superconducting systems, such as high-Tc cuprates (HTSC) 
[6-13], fuUerides [14,15], and two-dimensional (2-d) organic BEDT systems [16]. Universal 
nearly linear correlations have been found between Tc and the muon spin relaxation rate 
a{T — > 0) oc 1/A^ oc ng/m* (superconducting carrier density / effective mass) in the un- 
derdoped region of many HTSC systems and in some other exotic superconductors [6,7] 
Such correlations are seen also in HTSC superconductors having extra perturbations, such 
as overdoping [8,12,17], (Cu,Zn) substitutions [18], or spontaneous formation of nano-scale 
regions with static stripe spin correlations [19,20]. In all these systems, Tc follows the cor- 
relations with superfiuid density found for less perturbed standard HTSC systems. These 
results indicate that the superfiuid density is a determining factor for Tc in the cuprates [21]. 

In general, a strong dependence of Tc on the carrier density is not expected in conven- 
tional BCS theory [22] where Tc is determined by the mediating boson (phonon) energy 
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scale and the density of states of carriers at the Fermi level which govern the charge-boson 
(electron-phonon) coupling. For 2-d metals, the density of states does not depend on the car- 
rier density in the simplest case of non-interacting fermion gas. Therefore, the BCS theory 
has a fundamental difficulty in explaining the observed correlations. In contrast, an explicit 
dependence of can be expected for the condensation temperature Tg in Bose-Einstein 
(BE) condensation of a simple Bose gas, as well as for the Kosterlitz-Thouless transition 
temperature Tkt for a 2-d superfluid [23]. In BE and KT transitions, the transition tem- 
perature is determined simply by the number density and the mass, since the condensation 
is decoupled from the formation of condensing bosons in these two cases. The universality 
of the Tc vs ng/m* relationship observed beyond the difference of systems, such as cuprates, 
fuUerides, orangics, etc., may be related to this feature. Pictures proposed to explain the 
correlations between and the superfluid density in the cuprates include crossover from 
Bose-Einstein to BCS condensation [24-28] and phase fluctuations [29] . In the present work, 
we extend our study to intercalated nitride-chloride systems, seeking further insights into 
such phenomenology. 

We have also determined the upper critical field Hc2 of nitride-chloride systems from 
magnetization and resistivity measurements. Although it is not easy to determine i?c2 in 
layered superconductors due to the strong superconducting fluctuations, Hao et al. [30] de- 
veloped an approach to overcome such difficulty using a model for reversible diamagnetic 
magnetization of type-II superconductors which have high k values. Here, k is the Ginzburg- 
Landau parameter defined as the ratio of the penetration depth A to the coherence length ^. 
In this model, one calculates the free energy including the supercurrent kinetic energy, the 
magnetic-field energy, as well as the kinetic-energy and condensation-energy terms arising 
from suppression of the order parameter in the vortex core. This method has been suc- 
cessfully applied to various high-Tc cuprate superconductors. We will apply this model to 
superconducting Hf(Zr)NCl in the temperature and field region where the effect of super- 
conducting fluctuations is negligible. 
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II. SYNTHESIS AND CHARACTERIZATION 



The parent compounds HfNCl and ZrNCl were synthesized by the reaction of Zr or Hf 
powder with vaporized NH4CI at 600 °C for 30 minutes in N2 flow [31]. The resulting powder 
was sealed in a quartz ampule and was purified by a chemical vapor transport method with 
temperature gradient [32]. For the purification, an end of the ampule with pre-reacted 
powder was kept at 800 °C and the other end, where purified powder is collected, at 900 °C 
for three days. No impurity phase was detected in the purified powder from x-ray diffraction. 

Since a Li-intercalated sample is sensitive to air, intercalation was performed in a glove 
box. We used three intercalation methods [33] to prepare a variety of samples: (i) Li- 
intercalated ZrNCl samples were prepared by soaking parent ZrNCl powder in three kinds 
of butyllithium solution. We used 2.0 M (mol/l) n-butyUithium solution in cyclohexane, 
1.3 M sec- butyllithium solution in cyclohexane, and 1.7M feri-butyllithium solution in pen- 
tane for this purpose. In this order, reducing ability becomes stronger and therefore higher 
concentration of Li atoms can be intercalated into the parent compound. We used solution 
which contained butyllithium corresponding to more than 2, 5, and 5 equivalents of ZrNCl 
for n-, sec-, and tert-butyllithium, to avoid the dilution of butyllithium in the intercalation 
process, (ii) We performed co-intercalation of Li and organic molecules into ZrNCl by soak- 
ing Li intercalated ZrNCl powder (prepared by the method (i) using n-Butyllithium solution) 
in enough THF or PC. (iii) Li- and THF-co-intercalated HfNCl samples were prepared by 
soaking HfNCl powder in various concentrations (2.5-100 mM) of hthium naphthalene (Li- 
naph) solution in THF. Since phase separation was often observed in the samples prepared 
by the method (iii), we selected single phase samples for the measurements in this study via 
characterization from Tc as well as c-axis lattice constant. The Li- and THF-co-intercalated 
HfNCl sample with the largest c-axis lattice constant according to the method (iii) was pre- 
pared at Hiroshima University, while all the other samples at Columbia University. Table 1 
shows a list of these samples. 

The chemical composition was determined by inductively coupled plasma atomic emission 
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spectroscopy (ICP-AEM) and CHN elemental analysis. The powder samples were pressed 
into pellets under uniaxial stress and sealed in cells made of Kapton film and epoxy glue for 
x-ray measurements. The x-ray rocking curve of (00/) peaks for cleaved surface (inside of a 
pellet) shows a peak with half width at half maximum of ~ 8 degrees, which indicates that 
the bulk of the samples have well-aligned preferred orientation and are suitable for studies of 
anisotropic properties. Essentially similar rocking curves were also observed for as-prepared 
surfaces of pellets. 

The superconducting transition temperature of the samples was determined from mea- 
surements of magnetic susceptibility x- Fig- 2, we show the results of x in intercalated 
ZrNCl specimens which were examined by //SR measurements. The values of x = M/H 
are of the order for ideal perfect diamagnetism — S/Stt for spherical samples and —l/An for 
long cylindrical samples with the field parallel to the cylinder axis. The determination of the 
absolute values of x, however, can be affected by such factors as non-spherical sample shape, 
sample morphology and residual field in a SQUID magnetometer. These factors might have 
caused deviation of some of the zero-field cooled shielding values of x from —l/iir. 

In Table 1, we summarize the composition, synthesis method, distance between adjacent 
stacking units (1/3 of the c-axis lattice constant Cq, see Fig. 1), and the superconduct- 
ing transition temperature Tc for these samples. The stacking unit distance is 9.4 A for 
the Li-intercalated samples prepared by the method (i) without co-intercalation of organic 
molecules. This value is almost the same as that for unintercalated parent compound (9.3 
A) without Li. With increasing Li concentration x from 0.17 to 0.6 (in the samples without 
co-intercalation of organic molecules), decreases from 14.2Ktoll.7K. Figure 3(a) shows 
the X dependence of Tc in Li3;ZrNCl and Li3;HfNCl samples, with/without co-intercalation, 
obtained in the present study. In Fig. 3 and Table 1, we find that: (1) Tc shows a slow reduc- 
tion with increasing x; (2) for close x values, Tc does not depend much on the stacking-unit 
distance Co/3 (see Fig. 3(b)). These results are quahtatively consistent with the reported 
results in the Li, K, and Na doped ZrNCl samples with/without co-intercalation of organic 
molecules [34]. Decrease of Tc with increasing charge doping is reminiscent of the case of 
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overdoped high-Tc cuprate superconductors. 

To the best of our knowledge, this is the first report of success in Li-THF intercalation 
into ZrNCl by sec- and tert-butyllithium. The stacking unit thickness is 13.3 or 18.7 A for 
the methods (ii) and (iii). As has been reported, Tc is almost unaffected by the expansion 
of the stacking-unit distance from ~ 9.4 A to ~ 13.3 A for the samples with Li content 
X ~ 0.17 (see Fig. 2 (b)). Systems based on HfNCl has Tc = 25.5 K, nearly a factor of two 
higher than that for intercalated ZrNCl. 

III. //SR: EXPERIMENTAL 

Our iiSR experiments were performed at TRIUMF, the Canadian National Accelerator 
Laboratory located in Vancouver, Canada, which provided a high intensity and polarized 
beam of positive muons. Each pressed pellet sample, with the c-axis aligned, was sealed in 
a sample cell which has a Kapton window and mounted in a He gas-flow cryostat with the 
c-axis parallel to the direction of muon beam. Transverse external field (TF) was apphed 
parallel to the beam direction, while muons are injected with their initial spin polariza- 
tion perpendicular to the field/beam direction. Low-momuntum (surface) muons with the 
incident momentum of 29.8 MeV/c were implanted in the pellet specimens. The average 
stopping depth, 100-200 mg/cm^, assured that the majority of muons are stopped within 
the specimen, after going through Kapton windows of the cryostat and the sample cell. 

Plastic sintillation counters were used to detect the arrival of a positive muon and its 
decay into a positron, and the decay-event histogram was obtained, as a function of muon 
residence time t which corresponds to the time difference of the muon arrival and positron 
decay signals. The time evolution of muon spin direction/polarization was obtained from 
the angular asymmetry of positron histograms, after correction for the exponential decay 
exp(— t/r^), where r^^ = 2.2 /is is the mean lifetime of a positive muon. Details of /^SR 
technique can be found, for example, in refs [35-37]. 
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The asymmetry time spectra A{t) were fit to a functional form; 

A{t) = A{0) exp(-(72iV2) X cos (cut + 0) , 

where A{0) is the initial decay asymmetry aXt — 0. The muon spin precesses at the frequen- 
cies (jj — jfj,Hf,xt, where 7^ is the gyromagnetic ratio of a muon (2n x 13.554 MHz/kG) and 
Hext denotes the transverse external magnetic field. As shown in Fig. 4, for an example of 
Lio.iyZrNCl, this oscillation exhibits faster damping in the superconducting state due to in- 
homogeneous distribution of internal magnetic fields in the flux vortex structure. In pressed 
pellet samples of random or oriented powder, this relaxation can usually be approximated 
by a Gaussian decay which defines the muon spin relaxation rate a. 

For systems except for ZrNCl-Li-THF, the //SR results were analysed by assuming a 
single-component signal, which shows a reasonble agreement to the data as in Fig. 4(b). 
In ZrNCl co- intercalated with Li and THF, the relative value of the shielding susceptibility 
was significantly lower than those of other samples, as shown in Fig. 2. Although it is not 
clear, this reduced susceptibility could possibly imply a finite fraction of superconducting 
volume. As a precautionary measure, by fitting selected low-temperature signals in field- 
cooled and zero-field cooled procedures to an asymmetry function having two-component 
signals, we estimated an upper-limit of the relaxation rate for the superconducting fraction. 
This upper-limit is shown (in Figs. 6 and 7) by the error-bar placed to the right side of the 
main point for a which was obtained for a single-component asymmetry function. 

IV. /xSR: SPECTRA AND RELAXATION RATE 

Figure 4 shows time spectra of muon decay asymmetry for a representative sample above 
and below Tc. In the normal state above Tc, the oscillation shows a small relaxation due 
to nuclear dipole fields. We denote this relaxation rate as (T„. Below Tc, the relaxation 
becomes faster due to additionaal field distribution from the flux vortex lattice. For each 
specimen, the zero-field /iSR spectra obtained above and well below Tc did not show any 
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difference. This assures that the temperature dependence of the relaxation rate observed in 
TF is due to superconductivity alone, and also imphes that there is no detectable effect of 
time-reversal symmetry breaking, contrary to the case of UPta [38] and Sr2Ru04 [39]. 

The effect of the superconducting vortex lattice can be obtained by subtracting this 
normal-state background (T„ from the observed relaxation rate aob- Since the nuclear dipolar 
broadening and superconducting broadening of the internal fields do not add coherently, here 
we adopt quadratic subtraction to obtain the relaxation rate a due to superconductivity as: 



Note that this procedure makes the error bar rather large around a = 0. 

Figure 5(a) shows the temperature dependence of the relaxation rate a for the samples co- 
intercalated with organic molecules having expanded interlayer distance. With decreasing 
temperature, the relaxation rate begins to increase below the superconducting transition 
temperature Tc. At the flux-pinning temperature Tp, the zero-field-cooling (ZFC) curve 
begins to deviate from the field-coohng (FC) curve. In the ZFC procedures, flux vortices 
are required to enter the specimen from its edge and move a large distance before reaching 
their equilibrium position. Below the pinning temperature Tp, this long-distance flux motion 
could be prevented by the flux pinning, resulting in an highly inequilibrium flux lattice and 
more inhomogeneous field distribution at muon sites. This behavior has been observed in 
earher //SR studies of HTSC [40], BEDT [16] and some other systems. In both systems 
shown in Fig. 5(a), we find that Tp is much lower than Tc, which is a characteristic feature 
for highly 2-d superconductors. The /xSR results of Tp for the present systems agree well 
with those from magnetization measurements discussed in Section VII. 

We performed FC measurements using a wide range of external transverse magnetic 
fields H^xt, and found no significant dependence of a on H^xt from 40 G to 1000 G, as 
shown in Fig. 5. In TF-//SR measurements in highly 2-d superconductors, such as Bi2212 




and 
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or (BEDT-TTF)2Cu(NCS)2, application of a high external magnetic field transforms 3-d 
flux vortex structure into 2-d pancake vortices, since higher field implies stronger coupling 
of flux vortices within a given plane and higher chance for the vortex location in each 
plane to be determined by random defect position on each plane [41]. The absence of held 
dependence in our measurements implies that corrections for the 2-d vortex effect is not 
necessary in the present study. This situation is expected for our c-axis aligned powder 
specimens. The relaxation rate cr{T) shows a tendency of saturation at low temepratures in 
all of the measured samples of nitride-chloride systems in the present study. This behavior 
is characteristic for s-wave superconductors. However, experiments using high-quality single 
crystals are required for a conclusive determination of the superconducting pairing symmetry. 
In the case of HTSC cuprates, d-wave pairing was established only after //SR results on high- 
quality crystals of YBa2Cu30j^ became available [42]. 

Figure 5(b) shows the temperature dependence of a for the samples without organic 
co-intercalant. The dependence of a on temperature T and held Hg^t of the fleld-coohng 
results was essentially similar to that for specimens with co-intercalation in Fig. 5(a). In 
these systems, we determined the pinning temperature Tp by magnetization measurements 
instead of by //SR, and show the results in section VII. 

V. /xSR: COMPARISON WITH OTHER SYSTEMS AND SUPERFLUID ENERGY 

SCALES 

The ^SR relaxation rate due to the penetration depth is related to the superconducting 
carrier density n^, effective mass m*, the coherence length ^ and the mean free path / as 

, n 47re^ Ug 1 
a oc X = — — X — X 



g2 ^* I 

The proportionality to ng/m* comes from the fact that this effect is caused by the super- 
conducting screening current, and consequently reflecting the current density in a similar 
way to the normal state conductivity of a metal which is proportional to the carrier den- 
sity divided by the effective mass. As will be shown later, ^ is estimated to be 80-90 A 
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in the present nitride-chloride systems. The mean free path cannot be determined at the 
moment, since a high-quahty single crystal is not yet available. In this situation, we proceed 
with the following arguments by assuming that the system falls in the clean limit << /). 
Clean limit has been confirmed in many other strongly type-II superconductors, such as the 
cuprates and BEDT systems. 

In highly anisotropic 2-d superconductors, the penetration depth measured with the 
external field parallel and perpenducular to the conducting plane could be very different. 
For the geomentry with Hg^t perpendicular to the conducting plane, related to the in-plane 
penetration depth Xab as in the present study, the superconducting screening current flows 
within the plane, resulting in the more effective partial screening of H^xt and the shorter A 
compared to the case with Hg^t parallel to the planes. In the present work, our specimen 
has a highly-oriented c-axis direction within -|-/- 8 degrees, and we regard our specimen 
as equivalent to single crystal specimens in terms of anisotropy. A theory /simulation work 
[43] shows that for un-oriented ceramic specimens of highly 2-d superconductors, value of a 
should be reduced by a factor of 1/1.4 from the value for single crystals observed with Hg^t 
applied perpenducular to the conducting planes. 

In Fig. 6, in a plot of versus the low-temperature relaxation rate a{T — > 0) oc ng/m*, 
we compare the results of the present nitride-chloride systems with those from cuprate and 
organic BEDT superconductors. The point for the BEDT system was obtained in //SR 
measurements using single crystal specimens [16]. The shaded area denoted as cuprates 
represents the universal linear correlations found for un-oriented ceramic specimens of un- 
derdoped YBCO systems: we multiplied the relaxation rate in these YBCO by a factor of 
1.4 to account for the difference between single crystal and un-oriented ceramic specimens. 
The data points lie in possibly two different groups having different slopes in the T^/ a rela- 
tion. The first group with a higher slope includes the present nitride-chloride systems with 
co-intercalation of organic molecules and the BEDT system, all of which having highly 2-d 
character as demonstrated by their depinning temperature Tp being nearly a 1/3 to 1/2 of 
Tc. The second group includes YBCO cuprates and nitride-chloride systems without organic 
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co-intercalation, all of which have more 3 dimensional (3-d) character in the flux pinning 
property with Tp closer to T^. The irreversibility and depinning behavior in Lij-ZrNCl with- 
out organic co-intercalation was studied not by /^SR but by magnetization measurements as 
described in section VII. 

The relaxation rate observed by /xSR is determined by the 3-d superfluid density rig/m*, 
as this is a phenomenon caused by the screening supercurrent density in bulk specimens. 
With the knowledge of interlayer spacing Cint, one can convert 3-d density ns/m* into 2-d 
density on each conducting plane as nsidj"^* — ns/m* x Cint- For systems having double- 
layer conducting planes, such as the present nitride chlorides or some family of the cuprates, 
the average interlayer spacing depends on whether or not the double layer is regarded as 
a single conducting unit or two. In our previous reports for cuprates [26], we treated the 
double layer as two single layers. We shall follow this approach here, and define the Cint to 
be a half of the stacking unit distance as cint — co/6. 

In Fig. 7, we show a plot of Tc versus the 2-d superefluid density ns2d/nn* represented 
by the value a{T — > 0) x Cinf We include a point obtained in c-axis oriented ceramic 
specimen of YBa2Cu307 [44] (without multiplying a factor 1.4 to a since this specimen 
had an almost perfect alignment of c-axis direction). We find that most of the data points 
share a unique slope in Fig. 7. This result suggests two features: (1) within the nitride- 
chloride systems, 2-d superfluid density ng^dlnn* is a determining factor for T^^ and (2) the 
2-d superfluid density may even be a fundamental determining factor for T^. among different 
superconducting systems. However, the second conclusion (2) must be taken with caution, 
because, this analysis depends on our treatment regarding single versus double layers, and 
also because recent data on T12201 [8]. and Bi2201 [45] cuprates, having very large interlayer 
distance Cj„t >12 A, show universal behavior with the results from YBCO (q^j ~ 6 A) only 
in a 3-d plot hke Fig. 6 but not in a 2-d plot hke Fig. 7 [26]. In contrast, the conclusion 
(1) is more robust, since all the nitride-chloride systems have double conducting layers, and 
since the predominant 2-d character is consistent with the absence of dependence of Tc on 
interlayer spacing in nitride-chlorides shown in Fig. 3(b). 
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Since the Fermi energy of a 2-d metal is proportional to the 2-d carrier density n2d divided 
by the in-plane effective mass m*, as Tp = {h^TT)n2d/iri* , the horizontal axis of Fig. 7 can 
be converted into an energy scale representing superconducting condensate. This conversion 
from penetration depth to the superfluid energy scale was first attempted by Uemura et al. 
[7] in 1991 and later followed by other researchers, including Emery and Kivelson [29]. 

In order to do such a conversion, one needs to obtain absolute values of the penetration 
depth A from the relaxation rate cr. The numerical factor in this cr to A conversion in cr oc X"^ 
depends on models used for analyses of relaxation function line shapes, fitting range of data 
analyses, treatment of single crystal versus ceramic samples, and some other factors. The 
Gaussian decay, which fits most of the data from ceramic samples quite well, is significantly 
different from the ideal field distribution P{H) expected for a perfect Abrikosov vortex 
lattice in triangular lattice. So, using a theoretical second moment for P{H) in Abrikosov 
lattice is not necessarily appropriate for data analyses in real experiments. After various 
simulations and consistency checks, we decided to adopt a factor which gives A = 2,700 A 
for cr = 1 iis~^ for a triangular lattice. Note that this conversion is for a standard triangular 
lattice, contrary to the statements of Tou and coUabrators [4,5] who have erroneously cited 
that we calculated A for a square vortex lattice. Then we can derive nsidj"^* from the 
observed values of cr and known values of Cint- In the horizontal axis of Fig. 7, we attach 
the 2-d Fermi temperature Tpid corresponding to the 2-d superfiuid density obtained in the 
above-mentioned procedure. 

A 2-d superfluid of bose gas, such as thin films of liquid He, undergoes superfluid to 
normal transition via a thermal excitation of unbound fiux vortices, as shown by Kosterlitz 
and Thouless (KT) [23]. For paired fermion systems composed of n fermions with mass m, 
forming a superfluid with boson density n/2 and mass 2m, the Kosterlitz-Thouless transition 
temprature T^t becomes 1/8 of the 2-d Fermi temperature Tpid of the corresponding fermion 
system. In the KT theory, the 2-d superfiuid density at the transition temperture Tkt should 
follow system-independent universal behavior: namely, (/i^vr) x ns2d/^* at T = Tkt equals 
TF2d/8- This universal relation was first confirmed by an experiment on He thin films [46]. 
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In systems close to ideal Bose-gas, the 2-d superfluid density shows almost no reduction 
between T — and T — Tkt [46]. Thus, in such a case, we would expect the points in Fig. 
7 (based on ns2d/'rn*{T = 0)) to lie on the Tkt line. In thin films of BCS superconductors, 
the superfluid density shows much reduction from the value of T = to T = Tkt, and 
the "KT jump of superfluid density" becomes invisible. This corresponds to the situation 
where the points in Fig. 7 lie far in the right side of the Tkt line. In Fig. 7, most of 
the points lie about a factor of 2 away from the Tkt line. The linear relation between Tc 
and ns2d/'m'*(T = 0) suggests relevance to the KT transition, as pointed out by Emery and 
Kivelson [29]. However, the deviation from the Tkt line implies serious difference from the 
ideal KT situation. 



VI. MAGNETIZATION MEASUREMENTS: EXPERIMENTAL 

Magnetization measurements were performed using a SQUID magnetometer (Quantum- 
Design) at Columbia. Aligned pressed samples were sealed in quartz ampules. The raw 
response curve was corrected by the subtraction of the quartz background curve measured in 
advance. In the normal state of the superconducting Hf (Zr)NCl samples, as well as the parent 
compounds, weak-ferromagnetic behavior is observed up to room temperature. This weak- 
ferromagnetic behavior changes by the intercalation. Therefore, we subtracted the weak- 
ferromagnetic contribution. We estimated this by extrapolating temperature dependence, 
assuming the Curie- Weiss law (M — C/{T — 9)) and fitting the normal state magnetization 
in the temperature range of 2.5Tc <T< 5Tc. In this temperature range, superconducting 
fluctuations can be neglected and the temperature dependence is slightly concave. We only 
used the data with the extrapolated weak- ferromagnetic contribution smaller than 10 % 
of the diamagnetic magnetization, to avoid an error from the assumption of the Curie- 
Weiss temperature dependence. The model developed by Hao et al. [30] was applied to 
the analysis of the reversible region. In this model, reduced (dimensionless) magnetization 
M' = M/V2Hc(T) and field H' = H/V2Hc{T) scales as a single function that contains the 
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Ginzburg-Landau parameter k as a unique parameter. In our analysis, we optimized Hc{T), 
in addition to k as a parameter independent of temperature. Resistivity measurements were 
performed using a well-aligned pressed sample with four electrodes that is sealed in a cell 
made of Kapton film and epoxy glue. 

VII. MAGNETIZATION MEASUREMENTS: SUPERCONDUCTING 

PROPERTIES 

Magnetization measurements were performed in Lio.iyZrNCl and Lio.isTHFo.osZrNCl 
with magnetic fields applied parallel to the c-axis. Figure 8 shows the results obtained 
after the corrections for the quartz ample background and for the weak-ferromagnetic con- 
tribution. We note that a crossing point exists in the M{T) curve under various magnetic 
fields for each system, which is characteristic of quasi- two-dimensional superconductors [47]. 
There are reversible temperature regions where ZFC and FC magnetization curves overlap 
each other. Below a certain temperature (the pinning temperature Tp), ZFC and FC magne- 
tization curves deviate. We notice that the reversible region is wider for Lio.isTHFo.osZrNCl. 
This result is consistent with a picture that, by the expansion of the interlayer distance, the 
interlayer coupling become weaker and the pinning of the vortices becomes less effective. 

In the data analyses in the reversible region of Lio.iyZrNCl and Lio.isTHFo.osZrNCl, we 
confined to the temperature region apart from Tc{H) in order to avoid ambiguity due to 
the superconducting fluctuations. As shown in Fig. 9, the data scale quite well to Hao's 
model in the whole reversible temperature range below Tc for the both systems. For all 
the data, M' « H' and hence the demagnetization factor can be ignored. This analysis 
yielded values of k, ranging between 50 and 80 (see Table 1), which indicates that these 
compounds are extreme type-II superconductors. In Fig. 10, we show the values of the upper 
critical field Hc2j/ciT) obtained down to T = 2 K in this process using Hao's model. The 
temperature dependence of Hc2,//c{T) fits well to an empirical formula ifc2,//c(0)[l — (T/Tc)^], 
as shown by the dashed hues in Fig. 10. We emphasize that the low temperature hmit value 
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Hc2,//c{T — > 0) can be obtained almost without any extrapolation using this formula: the 
resulting values are shown in Table 1. The i?c2//c value of ~ 4-5 T in ZrNCl-Li-THF system 
is about a factor of 2 smaller than Hc2 ~ 10 T in HfNCl-Li-THF system reported by Tou 
et al [5]. These results might indicate that Hc2 roughly scales with Tc- A similar nearly 
linear relation between i?c2 and T^, can be found in the Hc2 values for high-Tc cuprate 
superconductors in the optimum doping region. 

The critical temperature Tc obtained using Hao's model is 14.9 K for both samples, which 
agrees with the estimate from the onset of diamagnetism due to superconductivity. We no- 
tice that at i7 = 55 kG above i7c2,//c(0), a diamagnetic behavior due to superconducting 
fluctuation was observed in magnetization as shown in Fig. 8. Similar results due to critical 
fluctuations have been reported in HTSC [48-50] and BEDT [51] systems. We obtained 
the in-plane coherence length ^abiS^) and the in-plane penetration depth Aa(,,M(0) using ex- 
pressions i?c2(0) = 0o/27r^(O)^ and k — X/^. These results are also summarized in Table 1. 
We note that ifc2,//c(0), and Cab{^) are almost unaffected by the interlayer distance. This 
agrees with the view that the essence of the superconductivity in Hf(Zr)NCl is dominated in 
Hf(Zr)-N honeycomb double layers. The values of the penetration depth determined both 
from //SR and magnetization show reasonable agreement, although the former is ~ 20 % 
smaller than the latter for both compounds. 

In magentization measurements (see Fig. 8) and /xSR measurements (see Fig. 5(a)), 
the results become history dependent below a pinning temperature Tp for a given exter- 
nal field Hg^f This feature can be expressed by defining the irreversibility field Hi^r for a 
given temperature T as H^xtiT — Tp) — Hi„. Figure 10 also includes Hi^r as a function of 
temperature, determined from magnetization and //SR measurements. The results obtained 
from the two different techniques exhibit excellent agreement. ZrNCl superconductors have 
a quite large reversible region in the H-T plane. The temperature dependence of the ir- 
reversibility field fits well to a functional form Hirr{T) — Hirr{0){Tc/T — 1)", obtained for 
three-dimensionally fluctuating vortices [52], with n = 1.5 at low fields below H ~ 0.4 Tesla. 
This provides support to our assumption of 3-d vortex lines which we adopted in our anal- 
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yses of TF-/iSR spectra taken below H — 0.1 Tesla. At higher fields, the fitting becomes 
worse, similarly to Ref . [52] . This may be related to a dimensional crossover from 3- to 2-d 
vortex fluctuations. More careful measurements are necessary to conclude this point. 

In order to provide a cross-check for the results of Hc2{T), we performed magnetoresis- 
tance measurements on Lio.iyZrNCl. The temperature dependence of resistivity for two sets 
of field and current configurations is shown in Fig. 11. High resistivity of the order of 100 
mQcm and negative slope of the resistivity in the normal state at low temperatures could be 
due to grain boundaries and may not be intrinsic. The observed superconducting transition 
is broadened by superconducting fluctuations, weakly superconducting regions such as grain 
boundaries, and by vortex motion due to Lorenz force. We notice that the resistive broaden- 
ing is slightly larger for H/ /c, which is a natural consequence of significant superconducting 
fluctuations only for H/ /c. We defined Tc{H) where resisitivity shows 50 % drop of the 
maximum value. Figure 12 shows Hc2{T) obtained in this procedure. These absolute values 
of i?c2(0) for H/ /c agree reasonably well with those from magnetization measurements, in 
spite of the unreliable definition due to the broad resistive transition. The difference between 
the temperature dependences of the magnetization (Fig. 10) and resistive (Fig. 12) Hc2 data 
may be due the above-mentioned limitations of the resistive measurements. The anisotropy 
ratio of the upper critical field i^c2,±c/-f^c2,//c is roughly 3 as shown in Fig. 12. Although we 
do not have data for the system with co-intercalation of THF or PC, the anisotropy ratio 
would presumably increase in more 2-d systems with larger stacking unit distance. 

VIII. DISCUSSIONS AND CONCLUSIONS 

The quasi-two-dimensional nature of the superconducting state appears in various super- 
conducting properties of intercalated Hf(Zr)NCl. Tc correlates with a 2-d superconducting 
carrier density ns2d devided by effective mass m* rather than the 3-d counterpart. Dia- 
magnetic magnetization due to superconducting fluctuation for H/ /c is observed at high 
temperatures and high fields. The crossing point exists in M{T) curves measured at various 
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fields. The reversible region of magnetization becomes larger with the increase of interlayer 
distance, suggesting weaker interlayer coupling. 

In addition to these results, we note that Tc, ns2d/^* i and ^ah [iidij jcj^ all show moderte 
dependence on chemical doping level which presumably represents the in-plane normal-state 
carrier concentration, while remaining almost independent of the stacking unit distance. 
These parameters are closely related to the superconductivity mechanism of this layered 
superconductor. Since the coherence length is a measure for the pair size, independence of 
^ah on interlayer distance implies that interlayer coupling does not affect the pair formation. 
It is then possible to consider a picture in which fluctuating superconductivity exists within 
a given layer, while the layers are coupled weakly by Josephson coupling to achieve 3-d bulk 
superconductivity. 

On the other hand, \ah and Tp are strongly affected by the interlayer distance. The reduc- 
tion of \ah with increasing stacking-unit distance can be understood as a simple reduction 
of the supercurrent density caused by lower density of the planes. The strong dependence 
of Tp on Cint is not surprising: this behavior has been seen in many HTSC cuprates. 

Our results show that Hf(Zr)NCl with variable interlayer distance as well as carrier 
concentration is suitable for systematic studies of layered superconductors. In addition, low 
Hci value of this compound makes it easier to cover the whole superconducting region in 
the H-T plane, and helps our study of vortex phase diagram. 

In Fig. 6 and Fig. 7, we have compared the results from the nitride-chlorides with 
other layered superconductors. All the arguments in the previous paragraphs, as well as 
Fig. 7, give an impression that 2-d properties are predominantly important factor of all 
of these layered superconductors. However, comparison among different cuprate systems 
having different interlayer spacing Ci^t indicates that the 3-d interlayer coupling plays a very 
important role in determing T^. in the cuprates. Furthermore, the observed results in Fig. 
7 show about a factor 2 deviation from the Tkt hne. These results indicate that a simple 
theory for KT transition, whose = T/^y is unrelated to the interlayer coupling, is not 
applicable either to the cuprates nor to the nitride-chloride systems. Further experimental 
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and theoretical studies are required to determine the origin of this deviation. Studies of 
crossover from Bose Einstein to BCS condensation, in the case of 2-d systems, might provide 
a clue for understanding this feature. 

The absolute values of the penetration depth A, obtained from the /xSR and magneti- 
zation measurements, show a reasonable agreement. We notice, however, about 20-30 % 
difference in the values from the two different methods (see Table 1). /xSR and magnetiza- 
tion estimates of A often exhibit some disagreement of this magnitude, as can be found also 
in the cases of HTSC and organic systems. Ambiguity of A with 20-30 % would, however, 
correspond to ~ 50 % ambiguity in the estimate of the superfluid density. In this situation, 
it would be ideal if there were a method to cross-check the superfluid density derived from 
^SR results in a completely different perspective. 

In Bi2212 cuprate systems, Corson et al. [53] measured the frequency dependent super- 
fluid response, and found temperature Tkt above which the superfluid density depends on 
the measuring frequency. The superfluid density observed at T = Tkt was consistent with 
the value expected in the universal argument of KT. This provides an excellent system- 
independent calibration to the superfluid density. The difference between the superfluid 
density at Tkt and at T — > should correspond to the distance (in the horizontal direction) 
of the corresponding data point in Fig. 7 from the Tkt line. The /xSR Bi2212 data point 
in a plot hke Fig. 7 lie about a factor 2-3 away from the Tkt hne. This factor agrees 
reasonably well with the reduction of the superfluid density from the T = value to the 
T = Tkt value observed by Corson et al. [53] in Bi2212 system in a similar doping region. 
This satiafactory cross-check for the Bi2212 system indicates that our choice of the conver- 
sion factor between a and A was reasonable, and the superfluid density derived by //SR is 
very reliable. Of course, comparisons among /xSR data for different systems in a relative 
scale can be performed without being affected by an ambiguity of their absolute values of 
the superfluid density. 

We performed fiSK measurements on three different specimens based on ZrNCl with Li 
concentraitons 0.15, 0.17 and 0.4, and found that the results of 2-d superfluid density ns2d/i^* 
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for these systems do not show much difference among one another. This phenomenon could 
be explained by two different possibilities: (a) not all the Li atoms donnate carriers on 
the ZrN planes, and the Li concentration x does not serve as an indicator of normal-state 
carrier concentraion; or (b) all the Li atoms donnate electrons to the ZrN planes, but only 
a finite fraction of those normal-state carriers participate in the superfiuid. The situation 
(b) is similar to the case of overdoped HTSC cuprates [21,54], where an energy-balance in 
the condensation process seems to determine the superfiuid density. Further experiments 
on normal-state transport properties are required to distinguish between (a) and (b) in the 
nitride-chloride systems. 

In conclusion, we have synthesized and characterized several different specimens of in- 
tercalated nitride-chloride superconductors, and performed /xSR and magnetization mea- 
surements. The superconducting transition temperature and the upper critical field 
Hc2,// exhibit a nearly linear relation with the 2-d superfiuid density ns2d/'fn* i while show- 
ing almost no dependence on the stacking unit distance. These features suggest a highly 
two-dimensional nature of superconductivity in the nitride-chloride system. 
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FIGURES 




(a) 



(b) 




FIG. 1. Schematic figures of the crystal structure of Hf(Zr)NCl (a) without intercalation and 
(b) co-intercalated with Li and THF. The stacking unit thickness (= co/3, where cq is c-axis lattice 
constant) is also shown. 
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FIG. 2. Magnetic susceptibility x vs. temperature measured in H = 5 Oe in the field-cooling 
(FC) and zero-field-cooling (ZFC) procedures in the specimens of intercalated ZrNCl which were 
used in the //SR measurements. 
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FIG. 3. Dependence of Tc, as determined from the susecptibility, on (a) Li concentration x, 
and on (b) stacking-unit-distance Co/3 in Li^ZrNCl and Li2,.HfNCl with/without co-intercalation 
of THF or PC. 
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FIG. 4. The time spectra of muon asymmetry A{t) measured in Lio.irZrNCl at (a) T = 19 K 
(above Tc) and (b) T = 2.5 K (well below Tc) under the transverse external field of 500 G. The 
apparent precession frequency in this graph is modified from the actual precession frequency by 
the use of a rotating reference frame. The solid lines show a fit to Gaussian decay envelope. 
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FIG. 5. Temperature dependence of the muon spin relaxation rate o"(T) in (a) Hf(Zr)NCl 
co-intercalated with Li atoms and THF, and (b) Li^^ZrNCl without co-intercalation. The upper 
panel (a) shows the results obtained in the field cooling (FC) and zero field cooling (ZFC) proce- 
dures, while the results in the lower panel (b) were obtained in the FC procedure. Dashed lines 
are guides to the eyes. 
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FIG. 6. Correlations between Tc and the muon spin relaxation rate a{T 0) of intercalated 
Hf(Zr)NCl (present work), high-T^ cuprates [6] and (BEDT-TTF)2-Cu(NCS)2 [16]. The horizontal 
axis is proportional to the 3-dimensional superfluid density n^/m* in the ground state. The results 
of a for the cuprates, obtained using un-oriented ceramic specimens, are multiplied by a factor ~ 
1.4 for the comparison with those from nitride-chloride and BEDT systems obtained using single 
crystals and oriented ceramic specimens. 



34 



100 




0(T^O)xCi„t (Aiasec ) 



200 



400 



600 

TF2d(K) 



1600 



2000 



2400 



FIG. 7. Correlations between Tc and a{T 0) x Cmt of intercalated Hf(Zr)NCl, organic BEDT 
[16] and YBaaCugOT (YBCO) [44], where Cint stands for average interlayer distance. We regard 
the double layers in the nitride-chlorides and cuprates as two separate layers, and thus assume 
Cint = Co/Q for the nitride chlorides and Cint ~ 6A for YBCO. The horizontal axis is proportional 
to the 2-d superfluid density ns2d/'^*- To the horizontal axis, we also attach the corresponding 
energy scale, 2-dimensional Fermi temperature Tp2d^ obtained from the 2-d superfluid density. 
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FIG. 8. Temperature dependence of the magnetization of (a) Lio.iyZrNCl and (b) 
Lio.isTHFo.osZrNCl, in external magnetic fields applied parallel to the c axis. For the plotted 
data, the quartz ampule background is corrected and the weak-ferromagnetic contribution is sub- 
tracted. 
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FIG. 9. Magnetization as a function of applied field (parallel to the c-axis) in (a) Lio.irZrNCl 
and (b) Lio.isTHFo.osZrNCl, shown in the reduced (dimensionless) units. The solid curves represent 
a fit to the model of Hao et al. [30] 
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FIG. 10. Temperature dependence of the upper critical field ^^c2,//c find the irreversibility field 
Hirr in Lio.irZrNCl and Lio.isTHFo.osZrNCl. H^2,//c was obtained from magnetization measure- 
ments (abbreviated as M in the figure), and Hirr from magnetization and TF-^uSR measurements. 
The broken lines show a fit to the functional form of Hc2 given in the text. 
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FIG. 11. Temperature dependence of the resistivity of Lio.irZrNCl under external magnetic 
field for two sets of the field and current directions, meaured using a four-probe contact. The 
applied fields range from to 9 T, in the intervals of 1 T. 



39 



10 



Lio.nZrNCl 



H,2 (H//C) 



H,2 (H//ab) 



2 



4 - 



2 - 



10 



12 



14 



Temperature (K) 



FIG. 12. Temperature dependence of the upper critical field Hc2 of Lig.irZrNCl, obtained from 
magnetoresistance. 
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TABLES 

TABLE I. Synthesis method (see the text for details), stacking unit distance (co/3), supercon- 
ducting transition temperature Tc estimated from magnetization, magnetic field penetration depth 
\ah,iJ,SR{T — >■ 0) at zero temperature limit estimated from TF-/iSR, Ginzburg-Landau parameter 
K, upper critical field at zero temperature iJc2,//c(0)) coherence length at zero temperature ^06 (0)) 
and magnetic penetration depth at zero temperature \ab,M{^) estimated from reversible magne- 
tization for intercalated HfNCl and ZrNCl systems reported in the present work. TF-/iSR and 
magnetization measurements were performed under the magnetic field parallel to the c-axis. 



Sample Synthesis method co/3 Tc Aab,/iSR(r ^ 0) ac i?c2,//c(0) Cab(O) Aab,M(0) 

A K A T A A 

Lio.iyZrNCl (i) ra-Butyllithium 9.4 14.2 3700 56 4.7 83 4700 

Lio.4ZrNCl (i) sec-Butyllithium 9.4 12.5 3900 

Lio.eZrNCl (i) tert-Butyllithium 9.4 11.7 

Lio.isTHFo.osZrNCl (ii) THE 13.3 14.4 5200 76 4.2 88 6700 

Lio.isPCo.isZrNCl (ii) PC 13.3 14.6 

Lio.24THFo.i4HfNCl (iii) 8 mM Li-Naph 13.3 25.5 

Lio.rTHFo.sHfNCl (iii) 100 mM Li-Naph 18.7 25.5 3900 
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